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Abstract 

In this paper we propose a lattice analog of phase-space fractional Liouville equation. 
The Liouville equation for phase-space lattice with long-range jumps of power-law types 
is suggested. We prove that the continuum limit transforms this lattice equation into 
Liouville equation with conjugate Riesz fractional derivatives of non-integer orders with 
respect to coordinates of continuum phase-space. An application of the fractional Liou¬ 
ville equation with these Riesz fractional derivatives to describe properties of plasma-like 
nonlocal media is considered. 

PACS: 05.20.-y; 45.10.Hj; 61.50.Ah 


1 Introduction 

A consistent formulation of the nonlocal statistical mechanics was first constructed by Vlasov 
in BE]. Fractional calculus [31 HI 0 |6l 13 El [9] has a lot of applications in physics [TOl [HI [I2l 
H3l 14] [151 IB] and it allows us to take into account fractional power-law nonlocality of continu¬ 
ously distributed systems. Using the fractional calculus, we can consider fractional differential 
equations for conservation of probability in generalized phase spaces. The Liouville equation 
with phase-space fractional derivatives allows us to formulate fractional statistical mechanics 
that describes systems with power-law nonlocality. The fractional statistical mechanics can be 
considered as a special form of the nonlocal statistical mechanics [T3] - It should be noted that 
the use of the fractional derivative of non-integer order is actually equivalent to using an infinite 
number of derivatives of integer orders, which can be arbitrarily large values (see Lemma 15.3 
in [3]). 

We note that fractional Liouville and Bogoliubov equations are discussed in [IT] HU 19, 120] 
[2H [221123]. There are two different types of approaches to space-fractional generalization of 
the Liouville equation: 

1) The first approach to generalization of the Liouville equation in the framework of frac¬ 
tional dynamics has been suggested in P21 CES HU m\. This generalization is derived from 
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the normalization condition with the fractional integration over phase-space coordinates. This 
fractional normalization condition is interpreted as an equation in fractional-dimensional phase- 
space or in phase-space with fractional powers of phase-space coordinates. We should note that 
the fractional Liouville equations, which are suggested in [17J HU HU, [20], do not contain frac¬ 
tional derivatives of non-integer orders. 

2) The second approach, which is based on the Liouville equation with fractional-order 
derivatives, has been suggested in [2111221 13] . where the Caputo type of fractional derivatives is 
used. To obtain this equation, we use the conservation of probability for a fractional differential 
volume element of the phase space. The suggested fractional Liouville equation is used to 
derive the fractional kinetics equations. The Bogoliubov hierarchy equations with fractional 
derivatives with respect to phase space coordinates are derived in [21], [22], The Vlasov equation 
and the Fokker-Planck equation with the Caputo fractional derivatives are also obtained from 
the fractional Liouville equations. In [23j the Liouville equation with the Caputo fractional 
derivatives is used to describe media with spatial dispersion of power-law type that is considered 
in [49] . We should note that this phase-space fractional Liouville equation are considered only 
for Caputo fractional derivatives. 

As it was shown in [31] EU [33] (see also [311 ESI ESI EZ] and [3U ESI HU HD E2] EU El]), 
the continuum equations with fractional derivatives can be directly connected to lattice models 
with long-range properties. Long-range interaction and properties are important for different 
problems in statistical mechanics [2111251126] . in kinetic theory and non-equilibrium statistical 
mechanics [27] [23], in the theory of non-equilibrium phase transitions [22] 130] . A connection 
between the dynamics of lattice system with long-range properties and the fractional continuum 
equations is proved by using the special transform operation [31, 32] S3] and it has been applied 
to different subjects [38] 39] EU EH E2| E3, S3])- 

In this paper, we propose Liouville equation for lattice phase-space. This Liouville equation 
describes fractional dynamics of a distribution function on unbounded homogeneous lattice with 
long-range jumps from one site to other sites. We prove that continuous limit of the suggested 
lattice Liouville equation gives the fractional Liouville equation for continuum phase-space. 
This fractional Liouville equation contains generalized conjugate Riesz derivatives of non-integer 
orders with respect to coordinates and momenta. As an example, we consider an application of 
the fractional Liouville equation with these Riesz derivatives to describe properties of plasma¬ 
like nonlocal media. 


2 Liouville equation for lattice with long-range jumps 

In statistical mechanics [33], the basic concept is the ensemble that is a set of classical sys¬ 
tems identical in nature, which are subjected to forces determined by identical laws m, but 
distributed in phase-space. We can state that the statistical ensemble is a set of independent 
systems identical in equations of motion, but differing in their initial conditions [33] ■ It is 
usually assumed that the systems are continuously distributed in the phase-space. In this pa- 
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per, we consider systems of particles that are distributed in lattice phase-space. Note that the 
statistical method should be used not only for systems with a very large number of particles. 
It is necessary to use this method in every case, even that of a single particle in the simplest 
possible conditions [46]. In the Liouvillian picture, we can describe the dynamics of a statis¬ 
tical ensemble in the phase-space points (sites) through which the ensemble points move. In 
this case, properties of the ensemble are assigned to points (sites) in phase-space at each given 
time, without an attempt to identify the individual system of the ensemble. In this paper, we 
consider a statistical ensemble on the lattice phase-space. 

Let us consider a system of classical identical N spinless point particles, which are char¬ 
acterized by the phase-space coordinates q S) k and p s ^ with s = 1, k = 1,2,3, As a 

lattice analog of continuum phase-space for this system we will use an unbounded lattice in 
6A r -dimensional Euclidean phase-space M 6Ar . The lattice is characterized by space periodicity. 
For unbounded phase-space lattice we define QN non-coplanar vectors ai, ..., that are the 
shortest vectors by which a lattice can be displaced such that it is brought back into itself. For 
simplification, we consider a lattice with mutually perpendicular lattice vectors a^. We choose 
the Cartesian coordinate system such that the directions of the axes coincide with the vector 
a j, such that = a^ /1 a^ | are the basis vectors of the Cartesian coordinate system in M 6JV . This 
simplification means that we consider 6 A r -dimensional analog of the primitive orthorhombic 
Bravais lattice. 

Choosing the coordinate origin at one of the sites, all phase-space lattice sites can be num¬ 
bered by the vectors n g = (ni, ..., n 3 jy) and n p = (n 3 jv+i, • • •, n 6N ), where rij (j — 1,, 6N ) 
are integer. The position of a lattice site is defined by the vector 

6 N 

x(n) = x(n 9 ) + x(n p ) = n^, (1) 

3 = 1 

where 

3 N 6 N 

x ( n <?) = x ( n p) = r L a r ( 2 ) 

j =1 J=3JV+1 

Considering the statistical ensemble on the phase-pace lattice, the W-particle system is 
represented by a point-particle, which is moved through the lattice sites. This lattice particle, 
which is lattice analog of the Wparticle system, will be called the lattice quasi-particle. We 
assume that it can be localized in the sites of the phase-space lattice, i.e. the positions of 
the quasi-particle coincide with the lattice sites. Then the vectors n g and n p can be used to 
describe the quasi-particle. We will consider long-range jumps of the quasi-particle on the 
lattice phase-space. 

The distribution function of lattice quasi-particle, which is the probability density to be in 
the lattice site, will be denoted by p(n 9 , n p , t) = p(ni, ..., n 6 n, t), where the vectors n g and n p 
define the site. We can use the notation pi(n g ,n p ,t) instead of p(n q ,n p ,t) to emphasize that 
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it is a lattice function. The function p(n q , n p , t) satisfies the conditions 

+OO +OO 

E - E p(n 1} ... ,n 6N ,t) = 1, p(n 1 ,...,n 6 N,t)>0 (t E R). (3) 

n i=—oo riQN=—co 

In the lattice, any particle of the iV-particle system is described by the six numbers 

n(s) = (n s , n s+ 1 , n<j + 2 , n^N+si ^ 3 JV+s+i, ^37 v+s+ 2), (s = 1,..., N). (4) 


Since we assume that the system consists of the identical N particles, then the function 
p(ni,... , rtQNi t) has the symmetry with respect to permutations n(s) < —* n(s'), where 

1 < s, s' < N. We also can define reduced s-particle distribution functions by the equation 

p s (n q ,n p ,t) = Y ■ ■ ■ Y ( 5 ) 

n(s+l) n(A r ) 

where p s (n g , n p , t) depends on the vectors n(l),.. . n(s) only. In the case s = N, we can see 
that p N (n q , n p , t) = p(n q , n p , t). 

Let us consider the equation for distribution function on unbounded homogeneous phase- 
space lattice in the form 


dp(n q ,n p ,t) 

dt 


3 N 

E : 

3 = 1 


L,q 




Vj( m q , m p , t) p{ m,, m p , t)) + 


6N 

+ B l, p 

j=3N+l 




m P , t) p(m g , 



J ( n 9 , n p , t), 


( 6 ) 


where the lattice function p(n q ,n p , t) is the probability density to find the quasi-particle at 
site (n 9 ,n p ) at time t. In equation (JBJ), the term J(n p ,n p ,t) describes an external source. The 


lattice operators ID) L q 


and D LjP 


Pi 

3 


are defined by 




L J 


+oo 


/(m g , m p , t) = — Y K °- : , ( n 3 - m 3 ) /( m 9> rrl P- t), 


a 3 


rrin =—oo 


(7) 


where we have D, 


for j = 3 N + 1,..., 6N. 


L,q 

The 


for j = 1,..., 3 N and 


L,p 


with /L = an 


particle long-range jumps with length rij 


kernel K~.(rij — nrij ) of the operator (J7J) describes the quasi- 
rrij on the phase-space lattice. If K~.(rij — rn 3 ) > 0 
then the kernel characterizes the jump to the site with rij from all other sites with rrij ^ rij. If 
K~.{rij — rrij ) < 0 then the kernel describes the reverse process. In equation (JUJ) , the parameters 
aj (and [3j ) are positive real numbers that characterize how quickly the intensity of the jumps 
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in the lattice decrease with increasing the jump length | rij —rrij\. These parameters also can be 
considered as degrees of the power-law spatial dispersion in the lattice [38, SOj. Equation ([6]) 
is the lattice fractional Liouville equation that describes long-range jumps on 6 N -dimensional 
phase-space lattice. 

Let us consider the kernel K~.(rij — rrij) in the form 



">/! r (atj + 2 

ctj~+2 


3 a.j + 4 7r 2 (rij — rrij) 2 

2’ 2 ’ 4 


aq > -2, (8) 


where \F 2 is the Gauss hypergeometric function jST] . Note that expression (JHJ) can be used 
not only for aq > 0, but also for some negative values of —2 < oq < 0. In the kernel notation, 
the minus is used to mark that it is the odd functions of integer variable n £ Z such that 
K~ ( rij — rrij) = —K~ (rrij — rij ) for all rij , rrij £ Z and j — 1,... 6 N. 

The discrete Fourier transform 


K-ikj) = 


+oo 


—ikiriA ts— 


7ln=—00 


oo 

/y, (//;) = -2i V A' 


sm(fcinj) 

m =l 


(9) 


of the kernels K (rij) defined by (JHJ) has the form 


K aj (kj) = i sgn (kj) \kj\ aj . (10) 

Note that the function (IB]) can obtained by the equation 

K~ (n) = — — f k aj sm(nk)dk (11) 

71 Jo 

that is the inverse relation to equation (j9]). For integer values of oq, we can get simpler equations 
instead of (JHJ). For example, we have 


A7(n) = 


- 1 ) ? 


n 


K 2 (n) = 


( — l) n 7T 2(1 — ( —1)”) 


(-l) n 7T 2 


n 


irn J 


K^n) = 


n 


^ (12) 


n° 


where (1 — (—l) n ) = 2 for odd n, and ((—l) n — 1) = 0 for even n. 
In the general case, we can consider the kernels 


K a< (kj) = i sgn (kj \kj\ a > + o(|%|“ J ), (% -> 0), 


(13) 


where the little-o notation o(\kj\ aj ) means the terms that include higher powers of |Aqj than 
\kj\ aj . The form (fl3j) means that we consider lattices with weak spatial dispersion [38]. If we 
use condition (!H instead of ( 1101) . then we can consider wider class of kernels to describe the 
long-range lattice jumps. As an example of the kernel with (1 131) . we give 


(_l)(n+i)/2 ( 2 [(n + l)/2] - n) T(aj + 1) 
2“j T((aq + n)j 2 + 1)T((q! j - - n)j 2 + 1) ’ 
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where the brackets [ ] mean the integral part, i.e., the floor function that maps a real number 
to the largest previous integer number. The expression (2 [{n + 1)/2] — n ) is equal to zero for 
even n = 2m, and it is equal to 1 for odd n = 2m — 1. Note that the kernel (fT4l) is real valued 
function since we have zero, when the expression (—l )( n+1 )/ 2 becomes a complex number. It is 
easy to see that we can use equation (fT4l) for all integer values neZ. 

For wide class of physical iV-particles systems, we can use 


Vj(m q , m p , t) = Vj (nip), Fj (m,, m p ,t) = Fj (m,, t ). 


(15) 


These conditions mean that Vj are components of the lattice analog of the particle velocity 
vector, and Fj are components of the lattice analog of the force vector that is independent of 
the particle momenta. In this case, the Liouville equation © can be rewritten in the form 


d„(n,,n r ,t) + ^ Vj(m 


dt 


j= 1 


■p) 


L,g 


a, 


p(m 9 , m p ,t) + 


6 N 


j=3N+l 


Pj 

.1 


p(m 9 , m p , t) = J(n q ,n p ,t). 


(16) 


In the general case, the usual Leibniz rule for the lattice fractional derivative does not hold 
if oij 7 ^ 1. This means that we have the inequality 



/(m,)ff(m ? )) ^ /( m g) 


L,g 


Otj 


g(m q )+g(m q 


L,q 


a i 


/( m . 


(17) 


This property is analogous to the characteristic property of derivatives of all fractional orders 
and all integer orders aq ^ 1 [50]. 


3 Fractional Liouville equation for phase-space contin¬ 
uum 

In this section, we use the methods suggested in [31] [32, 33j [13] to derive a fractional Liouville 
equation for phase-space continuum with power-law non-localities. 

Let us define a transform operation that allows us to derive the fractional Liouville equation 
for continuum phase-space from the Liouville equation for lattice phase-space (l 6 l) . 

Definition The lattice-continuum transform operation Tl~>c is the combination 

Tl^c = ° Lim o Fa (18) 


of the following three operations: 
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1) The operation J 1 a is the discrete Fourier transform /^( n q ,n p ) —> d-A{/z/ng, n p )} = 
/(kg, k p ) of the lattice function /(rig, n p ) that is defined by 

+oo +oo 

/( k g, k p) = ^A{/L(ng,n p )} = ^ ... / L ((ng,n p )) e"* (k ’ x(n)) , (19) 

ni=—oo riQN =—oo 

where (k,x(n)) = (k g ,x(n 9 )) + (k p ,x(n p )), 

67V 37V 67V 

x ( n ) = 5^ n 3 a 3> x ( n ?) = n i a T x ( n p) = 

3=1 3=1 j=37V+l 

and aj = 2 tt/ kj 0 is distance between lattice particle in the direction a j. Equation ( T7P) ) means 
that we consider the lattice function /z,(n 9 , n p ) as discrete Fourier coefficients of some function 
/(kg, k p ) for kj e [-kjo/2, k j0 / 2], where / = 61V. 


2) The operation Lim is the passage to the limit /(k g ,k p ) —>■ Lim{/(kg, k p )} = /(k g ,k p ), 
where we use aj 0 (or kj 0 oo). It allows us to derive the function /(k g ,k p ) from 
f(K kp), where /( k g ,k p ) is the Fourier integral transform of the continuum function fc(c p p) , 
and /(kg,k p ) is the discrete Fourier transform of the lattice function /z,(n g,n p ) 7 where 


}l{ n g ,rip) 


' 67V 


fJ( 27 r/A: j o) 


u=i 


/c( x (rig), x(rip)), 


and UjOj = 2-KHj/kjQ —> qj for j — 1,, 3N and n 3 a 3 = 2nrij/kjo pj for j = 31V+1,. .., 61V. 

3) The operation T~ x is the inverse integral Fourier transform /(kg, k p ) —» d r_1 {/(kg, k p )} = 
/c/q, p) that is defined by 

1 /*+oo /»+00 

/c(q,p)=^ 1 {/(k a ,k p )} = — w / dk ,.../ dk m f( k„,k„). (20) 

d —oo a-oo 


For simplification, we will use notations /(n, 


q? Up) for /i(n q , n p ) and /(q,p) for /c(q,p). 
From the context it will be clear which function is considered. 

We use the lattice-continuum transform operation 7l->c f° r the lattice functions /(nig, nip, £), 
for the product of two lattice functions (/(nig, m p , t) p(m 9 , m p , t)) and for lattice fractional 

derivatives of the functions, i.e. g(n q ,n p ,t) = D p a f f(m q ,m pj t). The operation Tl~+c can 
be applied not only for lattice functions but also for lattice operators. The operation 72_>c al¬ 


lows us to map of lattice derivative 


into a phase-space continuum fractional derivative 
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that is defined in Appendix. 



Proposition The lattice-continuum transform operation Tl~>c maps the fractional Liouville 
equation for phase-space lattice fO) with the lattice operators & into the fractional Liouville 
equation for phase-space continuum 


dp(%P,t) y 

dt ^ C ' q 

3 =1 




^(q,P,t)p(q,P,t) + 


67V 

E “a, 

j=37V+l 


Pj 

j 


Fj(q,p,t) p{q,p,t)) = J (q, p, t), 


( 21 ) 


where 


"C,q 


and D 


C,p 


are the continuum fractional derivatives with respect to phase- 

space coordinate qj and pj that are defined by for 0 < ay < 1, by equation l83\) for aj > 1, 
and by equation UfR) for integer odd otj (see Appendix). The functions p(q, p, t), Vj(q, p, t), 
F(q, p,t) ; J(q, p, t) are defined by the equations 


Ptq.P.*) =TL^c(p(n„n r ,t)), J (q, p, t) = Tl^c (J(n„ n„, t )), 

Vj{q,p,t) = T L ^ c (Vj(n, I ,n p ,t)), F/(q,p,t) = T L ^c (Fj(n„ n r , ()). 


( 22 ) 

(23) 


Proof. 

Applying the discrete Fonrier transform J -a to the first term of the lattice Liouville equation 
), we obtain 

'dp(n q ,n p ,t)\ _ dJ r A (p(n q ,n p ,t)) _ dp(k q ,k p ,t) 


F A 


dt 


dt 


dt 


where 


p(kq, k p , t) = F a (p(n q , n p , t)). 

Analogously, we can see that 

'dp(n q , n p , t) \ _ dF L ^c (p(n,, n p ,t)) _ <9p(q,p,t) 


Tt 


L^C 


dt 


dt 


dt 


where p(q, p, t) is defied by ()22]h 

The discrete Fourier transform T A of the second term of (J6]), gives 


37V 


•^a(E d i 

o=i 


a,- 


V)!m, r m r ,t)p(m q , m,„<)) = 


(24) 

(25) 

(26) 





















3N I +oo 

Y1 e ~ ikj nj aj b l 

3 =1 \n j =-oo 


a.j 


kj( m g, m p, t)p(m q , m p , t) j = 


+00 +oo 


3 TV . 

-\ a i \ 

j=l J \rij=—oo rrij =—oo 


Y1 e ^ Uj ^ K <*i ( n i “ m i) ( V A m 'D ri V t)p( m g , nip, t)) I = 


3 TV 


+00 


+oo 


XI e lkjn ' jajK aj { n 'i) Y1 (Vk(m q ,m p ,t)p(m q ,m p ,t) 


■ ! a 7 

J=1 3 


e ~ikjmjdj | _ 


i n ■=— oo 


nij=—oo 


3N 


■ i a n 

3 = 1 3 


'y ^ aj (K aj (kjaj) ( ^(kg, k p , t) * p(kg, k p , t) 


(27) 


where K a ,{kj dj) is defined by (JHJ), and n'- = rij—rrij. The symbol ( * ) q denotes the convolution 
with respect to k q that is defined by the equation 

\ r+ko/2 3N 

f (k«j, k p , t) * g(k qi k p , t )) =/ JJ dk] f(k! k p , t) * g(k q - k' k p , t). ( 28 ) 

7 9 d-fco/2 i=1 

Similarly, the transform of the third term of (J6j) gives 

\ 6JV i , x 

p( m„ m p , t) ) = -j- ( Kp.{kj dj) (Vj(k q , k pj t) * p{ k q , k pj t)) J , 

/ j=3V+l a j ' V ' 


6 N 

J-a( £ B Z 
V?'=3iV+l 


ft; 


3 


where ( * ) p denotes the convolution with respect to k p that is defined by the equation 

-+ fc 0 /2 


(29) 


f(kq,k p ,t) * g(k q ,kp,t)) = TT dk'j f(k q ,k' p ,t) * g{k q ,k p - k' t). (30) 

J V J—kn /2 


P J-ko/2 j =3N+ 1 
As a result, the Liouvillc equation has the form 


dp(k q , k p , t) 
dt 


3 N 


k P , *) * d(k g , k p , t)J + 
3 =i 9 


6V 


(31) 


+ ^ ^ ^Pj ®j) (kg; kp, t) * p(kg, kp, f) J ./(k<j, kp, t), 

j=3AT+1 P 

where the symbols ( * ) q and ( * ) p denote the convolution with respect to k f; and k p respectively. 
Then we use 


K a . ( a.j kj) = i sgn (kj) \dj kj\ aj , (j = 1,..., 3N), 


(32) 
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(33) 


K o. (o,j kj ) = i sgn(fcj) \aj kj\^ j , (j = 3N + 1,... , 6N). 


The limit a,- —> 0 gives 


K aA k j) = lim —a~ K (kj a-j) i kj k, aj (j = l,...,3iV), 


“j ^-0 dj 


K. ikj) = lim -j- K , (kjaj) = ik s \kjf‘ 1 (j = 3N + 1,..., 6JV). 


■" , j 4 


(34) 


(35) 


The ”hat”-kernel K a .(kj) is the discrete Fourier transform of the kernel of lattice oper¬ 
ator. The equation that defines K~.(kj ) has the form 


T. A 


an 


/(m 9 , m p , f) = -a-K (kj a/) f(k q , k p , t ), 


(36) 


where /(k„ k p , t) = ^{/(m,, m p , t)}. 

The ” tilde”-kernel K~.(kj ) is the Fourier integral transform J 7 of the continuum derivative. 
The equation that defines K~ ( kj ) is 


J 7 





(kj dj) f(k qi kp, f), 
a i 


where f(k q , k p , t) = 7{f( q, p, f)}. 

As a result, the limit a,j —> 0 for (1BTT) gives the Liouville equation in the form 


(37) 


~ / 3N 

dp(k g' t k, ’’ t] + y, K-(kj) (n(k a , k p ,t) » p(k„ k„,«)) + 

3 = 1 
6iV 

+ (-^?( k <n k p^) *P( k <n k P ^)) = (38) 

j=3V+l P 

where 

*«,(%) =!%|%r>- 1 , (j = l,...,3iV), (39) 

Kj i (k j )=ik j \k j fi-\ (j = 3JV + 1.6JV), (40) 

and 

p(k g ,kp,f) = Limp(k 9 ,k p ,f) J(k 9 , k p , f) = Lim J(k 9 , k p , t). (41) 

The inverse Fourier transform of (1381) gives the fractional Liouville equation f[ 2 Tj) . As a result, 
we prove that the lattice fractional Liouville equation (JHD is transformed by the operation 7 l->c 
into the continuum fractional Liouville equation (1211) . 

This ends the proof. 
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4 Fractional Liouville equation for Hamiltonian systems 

For wide class of physical iV-particle systems we can use 

Vj(q,p,t) = y(p), (j = l,...,3 N) (42) 

Fj(q,p,t) = Fj(q,t), (j = 3N + 1,..., 6N), (43) 

For Hamiltonian systems with potential forces, we have 


V i(P) = ^-> = 3 N) 

I'jim - 0 = 31V + 1,...,6JV), 


(44) 


(45) 


where the Hamiltonian is H( q, p) = T(p) + [/(q), the function U( q) is the generalized potential 
of the force, and T(p) is the generalized kinetic energy term. For simple case, we have V)( p) = 
Pj/M. Conditions (1421) and (I43T) mean that V) are components of the particle velocity, and Fj 
are components of the force vector that is independent of the particle momenta. In the case 
fljgp and (1431) . the fractional Liouville equation for phase-space continuum can be written in 
the form 


dp{ q, P, t) 
dt 


3 N 


j= 1 


J p(q,p,f)+ 


6N 


J 


p(q,p,t) = J (q, p, t). 


+ 4) D c, P 

j=3N+l 

As a special case, we also can consider fractional Hamiltonian systems 


Pi 


Vj{ q, p,t) = 


"c,p 


Fj(q,P ,t) = -P Ci9 


Pj+3N 
j + 3N 

Otj-3N 


H a,p{ q,p), (j = 1,...,3N), 


’ q [j - 3N 

where H a ^( q, p) is the generalized Hamiltonian function 
Liouville equation can be represented in the form 

<9p(q,p,t) 


^a,p(q,p), (j = 3N +1,...,6N), 


(46) 
where 

(47) 

(48) 


In this case, the fractional 


dt 


+ {-ff».<j(q,p),p(q,p, *)}<*,<s = ^(q.p.i), 


(49) 


where { , } aj p is the fractional Poisson brackets 

3 N 

{/(q,p),^(q ! p)}a > p = 


3 = 1 


@3 N+j 

3N + j 


/(q,p4) 


V c,q 




s(q,p4)- 


li 


















y C,q 


a,- 




" c,p 


@3N+j 

3iV + j 


g(<hp,t) 


(50) 


If all atj = 1 and f3j = 1, then (l50j) is the usual Poisson brackets, and (1491) is the usual Liouville 
equation for classical Hamiltonian systems. 

If we use (n q , n p ) instead of (q, p), and Bjj g , B 2 p instead of in equations 

then we get the correspondent lattice analogs of the considered equations. 


5 Fractional Liouville equation for nonlocal media 

As an example of application of the fractional Liouville equation with conjugate Riesz fractional 
derivatives, we consider description of non-local plasma-like continuum. 

For simplification, we consider an isotropic collisionless nonlocal plasma-like media, where 
all aj = a and (3j = 1, J(q, p, t) — 0 and Vj( p) = Pj/M. In this case, the V-particle distribution 
function p(q, p, t) is the product of one-particle reduced distribution function pi(x, p, £), 


N 


p( q, p,t) = Rpi(x„p„«), 


(51) 


S =1 


where we use usual space coordinates x s = L)y=i e j x sj and p s = e jPsj instead of the 

generalized phase-space coordinates q = (qi, , q 3 iy) and p = (pi,... ,P 3 n)- Here e 3 are the 
basis vectors of the Cartesian coordinate system in M 3 . The distribution function pi(x, p, t) 
describes a probability density to find the particle in the phase volume d 3 xd 3 p. 

In this case, the fractional Liouville equation 
Pi takes the form 

3 


t + 5>,(p> 


3 = 1 


where we use 


"c, x 


v c, P 


P i 


for the one-particle distribution function 

(52) 


Y.F 3+1 (x ] t)p-= 0 , 

_d_ 
dp 3 ' 


_3N + j 


Let us apply this Liouville equation with space-fractional derivatives B Cj , 
erties of nonlocal media. 

In the absence of the force held, the Liouville equation (j52j) gives 

3 


(53) 

to describe prop- 


dp 

dt 


u C,x 


a 


J 


P = o. 


(54) 


+ Y1 ^(p) 

3 = 1 

The solution of this equation will be denoted by p 0 ( x )P, t). For a weak force held, we can use 
the charge distribution function in the form 


Pi — do + <5p ; 


(55) 
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where po is the stationary isotropic homogeneous distribution function unperturbed by the 
fields, and 8p describes the change of po by the fields. In the linear approximation with respect 
to held perturbation, we have 


88 p 

dt 


3 = 1 


a 


Sp + F 3+i (x,t)^)= 0. 


(56) 


If we consider plasma-like media, then the force is the Lorentz force F = qE + g[V, B], where 
q is the charge of a particle that moves with velocity V = p/M in the presence of an electric 
held E and a magnetic held B. 

In an isotropic media, the distribution function depends only on the momentum, i.e. po = 
Po(|p|). In this case, the direction of the vector “ p 0 coincides with the vector p = MV, 

and its scalar product with [V, B] is equal to zero. Therefore, the magnetic held does not affect 
the distribution function in the linear approximation (J56]). As a result, we have 


88 p 

~dt 


+J2 (^-(p) 

3 = 1 


" C,x 


a 


J 


8p + q Ej(x, t) 


dpo 

dPj 


= 0 . 


The Fourier transform with respect to space and time gives 

3 3 


i (5^ s § n (%) I k j\ ay 3 - u)Sp + q ( E . 

3=1 3=1 ' 


dpo 

dPj 


= 0 . 


If we take the X-axis along k, then we have k = (k x , 0, 0), and 


(57) 


(58) 


T.sga(k t ) \kj\ a Vj = 

3=1 

where k x > 0 and sgn(k x ) = 1, since we use e x = k/|k|, and k x = |k|. In this case, equation 
gives 


i(\k x \ a V x ^uj)8p + qJ2( E ; 

3=1 


dpo 

dpj 


= 0 . 


Then we have 


8p = 


i I \k x \ a V x - c o 


U E 

3=1 


dpo 
J dpj 


(59) 


(60) 


In an unperturbed plasma-like media, the charge density is equal zero. The charge density 
that is perturbed by the held is 


Pcharge Q 



-Too 


<5pd 3 p = iq 2 



+oo 


1 V (E ■ — 

\k x \ a V x -u V 3 dp j 


d 3 P, 


(61) 
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where p c harge is the bound charge density. The electric polarization vector P is defined by the 
relations 

div P Pchargei (62) 

which has the Fourier transform in the form 

^(k, P) Pcharge• (63) 


The polarization P defines the electric displacement field D by the equation D = £oE + P. 
where e$ is the electric permittivity. Let the field E be parallel to k. Then P be parallel to k, 
and 

P= (e|l(|k|)-e 0 )E, (64) 

where £||(|k|) is the longitudinal permittivity. Substitution of (t6T]) and (1641) into (|63|) gives 


(£||(|k|) -£ 0 ) (k, E) = —q 2 




d 3 p, 


(65) 


where (k, E) is the scalar product of vectors k and E. Since we take the A"-axis along the 
vector k, then E = ( E x , 0, 0), and 


(k. E) = k x E x = \k\E x , 


3 = 1 


where we can use |k| instead of k x . 
Using the distribution function 


Po{Px ) = 


+oo 


Po{\p\)dp y dp z , 


we get the following equation to calculate the longitudinal permittivity 

q 2 f + °° 1 dp 0 (p x ) 


£ ||(l k |) = £o 


\k x \ a p x /M -uj dp x 


dpxi 


( 66 ) 


(67) 


( 68 ) 


where \k x \ = |k| can be used. For isotropic homogeneous case, we can use an equilibrium 
distribution po(p x ). 

Let us consider a plasma-like medium with the equilibrium Maxwell’s distribution 


p ° (fe) = TSPST exp ’ 

where ks is the Boltzmann constant. Then 

dpo(Px) = 2p x N q / Pi \ 

dp x 0F(2Mfc B T) 3 /2 exp ^ 2 Mk B T) ’ 


(69) 


(70) 
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where N q is the total number of particles per unit volume. 
Using k x = |k|, equation (1681) can be rewritten in the form 


a 2 N 

£i|(|k|) = £ 0 + Q q 


2 M 


^+oo 


Px 


|k| 1+ “ - v / / 7 t (2 Mk B T) 3 / 2 J-ac p x — Moo/\k\ a — iO 
Using new variables 

PX £ _ 


— exp 


z = 


y/2Mk B T 

equation (1711) can be represented in the form 

£ ||(l k |) = £q + 


M oo 


2 k B T |k|“’ 


^+°o 


Px 


2mkBT 


dp x • (71) 


(72) 


|k| 1+ “ ^/nk B T ./-co z-(-i 0 


— dz, 


(73) 


where 



--- -dz = y/^+P.V. 

d) 



fe 




dz + 27r£e 


-U 


(74) 


It should be noted that k, x, k x and x 3 are dimensionless variables. 

We consider equation (173]) for two cases that are characterized by the large and small values 
of the variable £. 

1) For small values £ <C 1, we have 


P.V. 


,+oo xe - z * 


dz = P.V. 


+°° £ e -(y+S ) 2 


dy = 


/ +°o . . 

e ~ y2 (y-'a - 2£ 2 - ( y - 1 + 2 y)e + (2 - (4/3 )?/ 2 )£ 4 + ..)dy = 

-oo x / 

= —2-v/vr£ 2 + a/vt£ 4 + • • • (£ « 1), (75) 

where y = z — £, and we take into account that the integrals of the odd terms in y are zero. 
Substitution of (1751) and (1741) into (1731) gives 


£||(l k l) 


£o + 


g 2 N q 

k B T Ikl^ 


q 2 N q Moo 2 | q 2 N q M 2 uo 4 
k 2 B T 2 |k | 3 “ +1 + 4fc|T 3 |k | 5 “ +1 + ' 


(76) 


In the case ( C 1, the imaginary part of the permittivity £ii(|k|) is relatively small but it is 
not exponentially small, because of the smallness of the phase volume, where the condition 
|k \ a p x /M — lo = 0 holds. 

2) For large values £ 1, we have 


P.V. 


f+°° xe~ z2 

7-00 


dz 



i -z/i 


dz 



dz = 


15 





















(£ »i), 


= ~\pK — 




2e 



(77) 


where we take into account that the integrals of the odd terms in z are zero also. Substituting 
(1771 ) and (1741) into (1731) . we get 


£ ||(|k|) = e 0 - 


q 2 N q 

Mu> 2 


k 


Q —1 


3 q 2 N q k B T 
M 2 oj a 


k \ 3a ~ 1 + ... . 


(78) 


For Maxwell’s distribution, an exponentially small part of the charged particles has the velocity 
V x = cu/|k| Vt = y/k^TJm, where Vp is the average velocity of charged particles. Therefore 
the imaginary part of the permittivity £||(|k|) is exponentially small. 

Equations (176]) and (178]) can be used to obtain the scalar potentials of electric held for the 
fractional nonlocal plasma-like media [171 39). 


6 Conclusion 

In this paper, the Liouvillc equation for unbounded homogeneous phase-space lattice with 
long-range jumps is suggested. Using the methods proposed in [31] 02, [33], we prove that 
the continuous limit transforms the suggested Liouville equation for lattice phase-space into 
the fractional nonlocal Liouville equation for continuum phase-space. This fractional Liouvillc 
equation contains the generalized conjugate Riesz derivatives on non-integer orders with respect 
to phase-space coordinates. As an example, we consider an application of the fractional Liouvillc 
equation with the Riesz derivatives of non-integer orders to describe properties of plasma-like 
nonlocal media. 

We assume that the lattice fractional Fokker-Planck equation, which is suggested in [33], can 
be obtained from the phase-space lattice Liouvillc equation that is proposed in this paper. The 
lattice Bogoliubov hierarchy equations can be easily derived from the suggested lattice Liouvillc 
equation by using the reduced distribution functions. The correspondent lattice hydrodynamic 
approximation can also be obtained from the lattice Bogoliubov hierarchy equations. It allows 
us to formulate fractional-nonlocal statistical mechanics |2U 13] on the lattice phase-space. The 
suggested lattice-continuum transform operator allows us to have a correspondence between 
lattice fractional statistical mechanics and fractional nonlocal statistical theory. 

The fractional Bogoliubov hierarchy equations with derivatives of non-integer orders with 
respect to phase-space coordinates can be obtained from the suggested fractional Liouvillc 
equation with Riesz fractional derivatives in the same way as was done for the Liouville equation 
with Caputo derivatives in [21, 22]- The Vlasov equation and the Fokker-Planck equation 
with the fractional derivatives of the Riesz type can be also obtained from fractional Liouvillc 
equations [33] - All these equations form a basis of fractional nonlocal statistical mechanics. 
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Appendix: Continuum fractional derivative of the Riesz 
type 


Let us consider the continuum fractional derivative D 


c 


property 


of the Riesz type that has the 


J 7 D, 


"c 


J /(x) ) ( k ) = i sgn (kj) |%r i (J r /)(k) (aj > 0). (79) 


For 0 < Q.j < 1 the operator D c a ? can be considered as the conjugate Riesz derivative [5] 

with respect to Xj. Therefore, the operator (179]) can be called a generalized conjugate derivative 
of the Riesz type. 

For 0 < otj < 1 the fractional operator d79|) can be defined by the equation 


v c 


a,- 


L J 


/(x) = 


_d_ 

8Xn 


- Zj) f(x + (Zj - Xj) ej) dZj, (0 < a, < 1 ), 


(80) 


where ej is the basis of the Cartesian coordinate system. The function R a . (x) is the Riesz 
kernel, that is defined by 


R a (x) = 


7i 1 {oi j )\x 


CJL4—1 


aj ^ 2 n + 1 , neN, 


~ 7 i 1 (aj)\x\ aj 1 In |x| aj = 2 n + 1, n6 N. 

The constant 71 (aj) has the form 

2 ^vr 1 / 2 r(a/ 2 )/r((l - otj)/ 2 ) aj ± 2 n + 1 , 

(_l)(i-°y)/2 2 <u-V/2 r(a/2) T(1 + [aj - 1]/2) aj = 2 n + 1, 


(81) 


71 («j) = 


(82) 


where iV eN and aj G M + . Note the distinction between the continuum fractional derivatives 


v c 


and the Riesz potential consists in the use of \kj\ aj instead of |k| aj . 
or aj > 1 the fractional operator d79|) can be defined by the equation 


O, 


c 


OLi 


/(x) = 


d 


di(m, aj — 1 ) dxj 


(A™f)(x) dzj, (1 < aj <m+ 1), (83) 


where (A™u)(x) is a finite difference of order m of a function /(x) with the vector step z j = 


^3 7 


pN 


for the point x = x j e j 7 The centered difference 


(A”/)(x,) = X>l)' 


ml 


n =0 


n\ [m — n 


/(x- (m /2 - n) Zjej). 


(84) 
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The constant di(m,aj ) is defined by 


_7r 3/2 A m (a i )_ 

2“jT( 1 + a/2)T((l + atj)/2) sin(7ra/2) ’ 


where 


[m/2] 

A m (a,) = 2 ]T(-1)-‘ 

s=0 


ml 


s\(m — s)! 


(m/2 — s) aj 


for the centered difference (I84p . The constant di(m,cuj ) is different from zero for all olj > 0 in 
the case of an even m and centered difference (A ™u) (see Theorem 26.1 in |3j). Note that the 
derivative (1831) does not depend on the choice of m > otj — 1. Therefore, we can always choose 
an even number m so that it is greater than a\ — 1, and then we can use the centered difference 
(IM|i for all positive real values of ctj. 

For integer odd values of aj, we have 


“c 


2m + 1 
J 


/(x) 


9 2m +7(x) 

dx 2 r +1 ’ 


(m G N). 


(85) 


means that the fractional derivatives D 


c 


of the odd orders aj are local 


Equation 

operators represented by the usual derivatives of integer orders. Note that the continuum 


derivative 


"c 


2m 

j 


, where m G N, cannot be considered as a usual local derivative d 2rn / dx 2m . 


of the order 2m with respect to Xj. 
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